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Abstract 

Combining recent perturbative analyses on the static QCD potential and the quark 
pole mass, we hud that, for the heavy quarkonium states cc, bb and ti, (1) ultra- 
soft (US) corrections in the binding energies are small, and (2) there is a stronger 
cancellation of IR contributions than what has been predicted by renormalon domi¬ 
nance hypothesis. By contrast, for a hypothetical heavy quarkonium system with a 
small number of active quark flavors (re; ~ 0), we observe evidence that renormalon 
dominance holds accurately and that non-negligible contributions from US correc¬ 
tions exist. In addition, we examine contributions of renormalons at u = — 1. 
As an important consequence, we improve on a previous prediction for possible 
achievable accuracy of top quark MS-mass measurement at a future linear collider 
and estimate that in principle 20-30 MeV accuracy is reachable. 


During the past few decades, there have been signihcant developments in the anal¬ 
ysis of heavy quarkonium systems using perturbative QCD. Developments in compu¬ 
tational technology greatly advanced our understanding on the nature of quark masses 
and interquark forces. We anticipate that eventually these developments will deepen our 
understanding on the structure of perturbative QCD in more general contexts. 

Recently an important step toward this direction has been achieved. A computation 
was completed of the four-loop relation between the quark pole mass and the mass in the 
modihed-minimal-subtraction scheme (MS mass) [1]. This result, when combined with 
other known results such as the three-loop correction (03) to the static QCD potential 
hQCD(’") [2], sets our analysis at a new stage, namely, at full next-to-next-to-next-to- 
leading order (NNNLO) in terms of short-distance quark masses. It realizes a cancellation 
of infra-red (IR) dynamics at this order. 

In this hrst analysis we report what can be learned by combining existing results. 
In particular we compare the results of [31 sm to make clearer the nature of the per¬ 
turbative series of the heavy quarkonium energies, concerning (1) corrections from the 
ultra-soft (US) energy scale and (2) the renormalon dominance hypothesis. In addition, 
we examine contributions of an ultra-violet (UV) renormalon at u = — 1 and discuss 
possible contributions of an IR renormalon at n = -|-1. 

Motivations for performing such an analysis can be stated as follows. A few years 
ago, a convincing evidence has been presented for the existence of IR renormalons in 
the perturbative series of the energy of a static color source, which has an IR structure 
common to the quark pole mass [5]. Hence, it is among general interests how accu¬ 
rately the renormalon dominance picture holds for the quark pole mass. Furthermore, 
contributions of US corrections to the quarkonium energy have collected attention since 
long time El Eli- Despite an original expectation of being dominating at IR, there 
have been evidences that US corrections are moderate in size from comparisons of the 
perturbative predictions with experimental data for the bottomonium spectrum illlQ], 
phenomenological potential models of heavy quarkonia m, and lattice computations of 
Uqcd(’") [I2]- However, extraction of an accurate size of the US corrections still remains 
a challenge [13] . 

Important applications of this type of analysis include precise determination of the 
masses of the heavy quarks c, b and t from the energy levels of the lowest-lying heavy 
quarkonium states [T3|. (For earlier works, see [15] and references therein.) In this 
paper we apply our new understanding to a study of the possible achievable accuracy 
of top quark mass measurement expected at a future linear collider. Today, a precise 
determination of the top quark mass is highly demanded, for a precision test of the 
standard model of particle physics (SM) [TB], and also since the top quark mass plays a 
crucial role in the vacuum stability of the SM at a very high energy scale HZI. Hence, 
progress in our understanding of the heavy quarkonium states may lead to an access to 
deep aspects of the SM. 

The pole-MS mass relation can be expressed in a series expansion in the strong 
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Table 1: Exact result of for 0 < n; < 6 in the full theory, with n/^ = 1 heavy quark and 
ni massless quarks. We use eq. (p|) obtained by a fit of the results in [T8l dU [T]. An error of 
±21.5 is assigned to each value. 
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Here, m = denotes the MS mass renormalized at the MS mass scale; as{fi) = 

af''^{jjL) represents the strong coupling constant in the MS scheme, where rii is the number 
of massless quark flavors {ni = 3, 4 and 5 for the charm, bottom and top quarks, 
respectively); the renormalization scale /x is set to m. In most part of this paper, we use 
the coupling constant of the theory with rii flavors only as the expansion parameter. The 
coefficients di can be obtained from the corresponding mass relations in the full theory 
(with Uh heavy quarks and ni light quarks), respectively, by rewriting them in terms of 
the conpling constant of the theory with ni light quarks only. 

Let us hrst summarize the results of the previous analyses, on which our analysis is 
based. Refs. [31 m estimated dg on the basis of different assumptions, prior to Ref. [I], 
which accomplished the exact compntation of d^: 

• Ref. [3] required stability of the perturbative prediction for 2mpoie ± Vqcd(?") at 
relatively large r. Essentially the only assumption made is that US corrections in 
hgcol’") do not deteriorate perturbative stability (which holds up to NNLO) at 
NNNLO. 


• Ref. [1] assumed renormalon dominance in rupoie and Vqcd(?") and estimated their 
contribntions from the latter. Contribntion of US corrections in Vqcd(?") was snb- 
tracted in this estimate^ 

• The exact values of d^ are obtained combining the results of direct perturbative 
computations in [HI [191 [I]- 

Only the values for ni = 3,4,5 are presented explicitly in the hnal form in [T] (for 
the full theory with Uh = 1). Since we need the values for other n^s in our analysis, we 
derive the exact result of ds given as a cubic polynomial of ni as 

Cfuii theory = -0.67814 nf + 43.396 nf - 745.85 ni + 3556.5 , (2) 

where an error of ±21.5 is assigned to its value for each ni. We determined the last two 
coefficients of eq. ([2D by a £t using the results of [HI [H] in addition to the result of [T]. 

* Since US corrections in Vgcoir) do not contribute to the renormalon at m = 1/2, this manipulation 
is justified within the renormalon dominance hypothesis. 
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Table 2: Summary table of relevant estimates and exact results of d^. The first line shows the 
estimates based on stability of the perturbative prediction for 2mpoie + Vqcd(?'); tbe second 
line shows the estimates based on renormalon dominance hypothesis; the third line shows the 
exact results (converted to the values in the n/ flavor theory). 

For the reader’s convenience, we list the exact result of in the full theory in Tab. [T] 
using this formula for 0 < n; < 6 . 

As already mentioned, we convert the above formula using the coupling of the theory 
with Til massless quarks only as the expansion parameter. This gives 

Cconverted = -0.67814 uf + 43.396 nf - 745.42 m + 3551.1, (3) 

with the same error ±21.5. In the rest of the analysis, we use this for various n^s. 

In Tab. [2] we summarize the two estimates and the exact result for 0 < n/ < 60 
The relative accuracies are compared visually in Fig. [T] Overall, we hnd a reasonable 
agreement of the previous estimates and the exact results, with respect to the assigned 
errors. The relative accuracies of the estimates are also fairly good, at order 10% level. 
These features provide certain justihcation to the used assumptions in these estimates. 

Furthermore, we can make a closer examination. In particular, the central (optimal) 
values of d^^ in the table and hgure carry important information on the respective 
assumptions. We should note that the errors of df^ are only systematic and have no 
statistical nature. Hence, by carefully contemplating on the origins of these systematic 
errors, we can extract the sizes and signs of the systematic effects. The agreement with 
respect to the systematic errors is a necessary condition for the validity of our analysis 
given below. 

In the cases ni = 3,4, 5, corresponding to cc, bb, tt quarkonium states, respectively, we 
see a good agreement of the estimates by [3] with the exact values, whereas the estimates 
by HI are slightly larger. On the other hand, for smaller n/ = 0,1, which correspond to 
hypothetical heavy quarkonium systems, the agreement between the estimates by jl] and 
the exact results is fairly good, whereas the estimate by [5] for n; = 0 is slightly smaller 
than the exact value. From these observations we derive the following interpretation: 

• For Hi = 3,4, 5, (i) US corrections in Vqcd(^") are small, and (ii) there is a stronger 
cancellation of IR contributions than what has been predicted by renormalon dom¬ 
inance hypothesis. 

• For n; = 0, (hi) renormalon dominance holds more accurately, and (iv) non- 
negligible contributions from US corrections exist. 

t Since we use the converted its values for m = 3,4,5 listed in this table are different from 

TABLE III of [T]. In this sense, the comparison in TABLE III of [T] is not consistent, since da’s in the 
different definitions are compared. Numerically the differences due to different definitions are small, 
nonetheless. 
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Figure 1: Comparison of (ds — for the (converted) exact value of and the 

two estimates, where denotes the central value of 

We explain the details in the following. 

The renormalon dominance hypothesis assumes that the expansion coefhcient of the 
perturbative series is dominated by a factorial (~ n\) growth [2U] . 

dn ~ const. X (2/3o)” ^ n > 1, (4) 

which stems from the singularity at u = 1/2 in the Borel transform of the perturbative 
series, [z/ = /9 i/(2/9q), and fdi denotes the {i + l)-loop coefficient of the beta function of 
Q!s.] Contributions from the analytic part at u = 1/2 are neglected. The comparison 
between and the central values of d3®*[4] shows that the renormalon dominance 

hypothesis works better for smaller ni. This suggests that the above factorial growth 
overwhelms contributions from the analytic part as /9o(> 0) becomes larger for smaller 
m. 

Another source of rii dependence of the renormalon dominance resides in the series 

[201 [22] 

~ l/(l/ — 1) _ Z/(l/ — l)(z/ — 2) 

Fn = ^^ -7—T—V—T-TT 7—T—77—T-TT7—T-d7 

n + u [n + ujyn + u — \) [n + v)[n + v — l)[n + v — 2) 



in eq. (33) of [1|. The factor multiplies the right-hand side of eq. (jl]), giving 1/n 
suppressed corrections, so that it shows how the expansion coefficient approaches the 
asymptotic form at large orders (n S> 1)II Fig. [2] plots the series (|5]) in our case n = 3 

^ Contributions from the analytic part at m = 1/2 are not included in the series Fn, as they are 
suppressed exponentially. 
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for different n^s. They exhibit the tendency that 1/n suppressed contributions become 
more important for larger ni, although the hrst term (=1) is by far dominating. Both 
of these ni dependences in analytic and 1/n suppressed contributions have been taken 
into account in the error estimates of [1]. The former error enters as scale dependences 
in the analysis of [1] and is the main source of errors. 



Figure 2: Each term of eq. (P for n = 3 and different n^’s: /i = ^ ci, /2 = C 2 , 

fs = (n+u)(n+l-i)~(n+u- 2 ) ^3- The first term (/o = 1) is omitted, since it is by far greater. 

One may wonder if the UV renormalon at m = — 1 contained in the pole mass gives 
a significant contribution to the perturbative series of the pole mass. Based on an 
analysis in the large-/9o approximation, we estimate that the contribution of the u = —1 
renormalon to is fairly small compared to the errors of listed in Tab. El This 

is consistent, since the analytic part at m = 1/2 contributes dominantly to these errors, 
and the u = —1 renormalon belongs to the analytic part. The analysis also suggests 
that the u = —1 renormalon contribution is not a dominant component of the analytic 
part. Another important feature is that, since the UV renormalon is Borel summable 
and gives a well-defined contribution, as long as we obtain a converging series of a 
physical observable (such as the heavy quarkonium energy level), the contribution of 
the u = —1 renormalon to the error estimate becomes small (arbitrarily small unlike 
IR renormalons). Indeed contribution to the error is minor at our present perturbative 
order. We give details of the analysis of the u = —1 UV renormalon in the Appendix. 

Similarly there may be effects by the m = 1 IR renormalon contained in the pole 
mass, whose properties are less known. Known properties are as follows [21]. (a) It is 
induced by the non-relativistic kinetic energy operator Zl^/(2m); (b) It is not forbidden 
by any symmetry, and parametrically it possibly induces an order uncertainty; 

(c) It does not appear in the large-/3o approximation. With this limited knowledge, it is 
not easy to estimate contribution of the u = 1 renormalon in the estimate of ds in [1|. 
In principle, this contribution is exponentially suppressed in the estimate of the m = 1/2 
renormalon in the pole mass and is encoded in the scale dependence in the error estimate 
of df [4]. 
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We turn to the estimates of hs by [3], which incorporate the fact that cancellation 
of IR dynamics occurs beyond the renormalon dominance hypothesis. It can be un¬ 
derstood using the potential-NRQCD effective held theory [8], in which interactions of 
a heavy quarkonium and IR degrees of freedom are systematically organized in mul¬ 
tipole expansion in r. The leading order interaction is given by the interaction of an 
IR gluon with the total color charge of the heavy quarkonium, which vanishes for a 
color-singlet system. The corresponding contribution to the binding energy is given by 
an r-independent IR part of 2mpoie + Vqcd(’") [23]. The cancellation between 2mpoie 
and Vqcd(?") is not restricted to the renormalon part, and the analytic part at u = 1/2 
contains such contributions. 

In this general framework, the lowest order non-canceled IR contribution to the energy 
is given by a double insertion of the dipole interaction between the color-electric held and 
heavy quarkonium, expressed in terms of a non-local gluon condensate of the form ~ (r- 
E°' f-E°'). It is dominated by contributions from the US energy scale, and the perturbative 
evaluation of this condensate at log ag) and 0{aj) has been incorporated in Vqcd(?") 
in the estimate of d^. In principle, the u = 1 renormalon in the pole mass (if it exists) 
can ahect 2mpoie + Uqcd(^")- However, the large mass limit m —)■ cxo is taken in the 
analysis of [3], so that the u = 1 renormalon (order Aq^p/m) is suppressed compared 
to the M = 3/2 renormalon (order Hence, the estimate of d3®*[3] should not be 

ahected by the u = 1 renormalon. Furthermore, in estimating the ehects of taking the 
large mass limit are small for rii = 3,4, 5, compared to the real c, b, t-quark mass cases, 
hence, our discussion is expected to be valid for these real heavy quarkonium systems. 
In the cases ni <2 and rii > 6, perturbative analysis makes sense only in a hypothetical 
static limit m — )■ cxo, and our discussion is conhned to this limit. 

In perturbative QCD, instability against scale variation in IR region is manifest for 
all the physical observables, reflecting the blow-up of the running coupling constant at 
IR. For a “good” observable, generally scale dependence decreases as the order of per¬ 
turbative expansion is raised. Empirically this happens not only in the ultra-violet (UV) 
direction but also stability extends to IR region as the perturbation order is increased. 
In the case of the heavy quarkonium energy, the leading source of IR instability is the 
non-local gluon condensate dominated by US corrections. The (optimal) values of the 
estimates of d^ in the hrst line of Tab. [2] are chosen to optimize the stability of the per¬ 
turbative prediction for the energy in the IR region at NNNLO. A very good coincidence 
of these values with the exact results for = 3,4 suggests that the US corrections are 
small for these systems. Here, we may set the criterion for “large” or “small” by whether 
the corrections deteriorate stability of the perturbative prediction or not. 

As shown in j3|, perturbative stability of 2mpoie + Vqcd(’") is sensitive to the precise 
value of ds, and this sensitivity turns out to be asymmetric with respect to the sign of 
a variation of dsJl If d^ is larger than a certain critical value, stability of the prediction 
is lost very quickly. This leads to a fairly sharp upper bound on the estimate of ds 
for each ni. By contrast, stability of the prediction is degraded only gradually if ds 
is lowered from its optimal value. In this regard, a marked result is that in the case 
ni = 0 the exact value of ds is on the verge of or slightly above the upper bound of ds 
required by stability of the energy. Since US corrections are expected to be the source 

^ Qualitatively the same feature is observed for the heavy quarkonium energy levels |24| . 
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of IR instability of the energy, we infer that the US corrections are sizable in this case. 
Oppositely, in the case ni = 5, the exact value of lies slightly below that required by 
optimal stability of the energy. Hence, in this case US corrections do not deteriorate 
perturbative stability in any essential way, and US corrections may well be regarded 
as “small.” Such a dependence of IR stability on ni may result from the fact that the 
running coupling constant blows up most rapidly for ni = 0, while the running becomes 
milder as rii increases. (Note that we consider rii massless quarks.) If an IR catastrophe 
of perturbative stability should ever occur, it would be expected to appear first in the 
most rapidly running case. To demonstrate explicitly the level of instability in the case 
ni = 0, we show a plot according to the analysis of [3]. Fig. |3]shows the scale dependences 
of 2mpoie + Vqcd(?") at a relatively large r, where perturbative stability up to NNLO is 
close to marginal. The NNNLO line is flatter than the NNLO line in the large fi region, 
however, it grows in the small fi region and starts to show a sign of instability. See [3] 
for more details of the analysis method. 



Figure 3: Scale dependences of 2mpoie + VQCD(r) at different orders of perturbative expansion 
in the case n/ = 0. Input parameters are 0^(3GeV) = 0.2 [as(/u) blows up at p 0.62 GeV], 
m = 16 GeV (a large value is chosen to suppress sub-leading renormalons in mpoie), and 
r = 0.5 GeV“^. The exact value of is used. A horizontal line is shown as a guide. 

There is a difficulty in quantifying the size of US corrections more directly. By 
definition, the US corrections are dependent on the factorization scale p/, which should 
satisfy the condition [8] 


, CaCHs 1 /„\ 

Aqcd, — < p/ < (6) 

2r r 

where Ca = Nc = 3 is the Casimir operator for the adjoint representation. Given the 
different rii dependences of and we conhrm that a simple logarithmic 

dependence of the US corrections on p/, proportional to Q!glog(p/r), cannot explain 
the difference, even if we assume a reasonable ni dependence of p/. This is expected, 
since except in dimensional regularization, which conceals power-like dependences on 
the scales, we expect a much stronger dependence ~ p^r^ of the US corrections. This 
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dependence should eventually turn into a dependence on the physical US scales, namely 
fXf should be replaced by Cac^s /and Aqcd, where presumably the latter is more 
dominant at larger This requires (at least) an analysis analogous to that of |1] 
incorporating the m = 3/2 singularity in addition. Furthermore, we would need to 
separate UV and IR contributions in perturbative expansion systematically, to be able 
to accurately extract the US contributions [2S1 [2S]- Such a detailed analysis is beyond 
the scope of this paper. 

Thus, for the case ni = 0, we are (for the time being) content with the observation that 
everything is consistent. The renormalon dominance hypothesis can accurately estimate 
ds by the method of [1]. As mentioned, it is plausible that the renormalon dominance 
works most accurately in this caseQ On the other hand, stability of 2mpoie + Vqcd(?") 
at IR can in principle be jeopardized by US corrections and, if at all, this is expected to 
happen for smaller n/. 

In contrast, for ni = 3, 4, 5, the analytic part at u = 1/2 has a larger relative signif¬ 
icance, and the central values of the estimates by jl] depart from the exact values; see 
Tab. [2] and Fig. [T] We can circumvent this problem in the method of j3], since cancel¬ 
lation of IR dynamics takes place in the analytic part as well, and the contribution of 
US corrections is expected to be milder than the ni = 0 case. Thus, we are led to the 
interpretation as presented in the beginning of this discussion. 

On the basis of our understanding up to this point, we reexamine the prediction of 
the energy level of the (would-be) toponium IS state, using the NNNLO formula for 
the IS* * energy level [271 [281 [29]. We compare with the analysis [30], which examined 
the IF energy level calculated in terms of frit = and in the £ expansion [3T] . 

The large-/3o approximation (a crude approximation based on renormalon dominance) 
was used for estimates of 03 and dsl/l We replace them by the exact values. The essence 
of the analysis [30] is to use the renormalon dominance hypothesis for estimating a 
perturbative error in the top quark mass determination from the energy level of the 
toponium IF state. As a result, about 40 MeV for an expected accuracy was predicted 
for determination of the top quark MS mass. 

All the qualitative argument of [30] based on renormalon dominance hypothesis 
should be valid, since, as we have verihed, the renormalon dominance is qualitatively 
a good approximation. Nevertheless, according to our above understanding, the accu¬ 
racy of the prediction is expected to improve, since the cancellation of IR dynamics 
occurs at a deeper level than that of the large-/9o approximation. In the tt system, the 
leading non-canceled IR contribution from US corrections is expected to be “small” if 
our understanding is consistent. 

^ /i/ dependence is canceled in physical observables. Hence, we are ultimately interested in the 
dependence of physical observables on the physical US scale. At lower orders of perturbative series, 
only the scale Cac^s /is visible. As the order is raised, perturbative expansion becomes more 
sensitive to the Aqcd scale. The leading dependence of 2mpoie + FQCD(r) on Aqcd should appear as 

j,2 

■ 

* This feature appears to be slightly reinforced for n/ = 0 by a cancellation of the contribution from 
the u = —1 UV renormalon and other contributions from the analytic part at u = 1/2; compare Tabs. [2] 
and[3l 

1 More accurately, a Fade estimate of 03 was used and the prediction of the 15 energy level was shown 
to be quite close to that of the large-/3o approximation. Since the difference is minor and irrelevant in 
our context, we refer only to the large-/3o approximation. 
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Figure 4: Comparison of scale dependence of the toponium IS energy at NNNLO from the 
previous analysis (the large-/3o approximation [30]) and that using the exact values of 03 and 
ds- The input values are mt = 165 GeV, as{Mz) = 0.1185, and ni = 5. Horizontal dashed 
lines are shown as a guide. 


Fig. m compares the scale dependence of the toponium IS energy by the previous 
analysis ED] and that using the exact values of 03 and d^. A marked difference is that 
the former prediction is much more unstable in the IR region than the latter. This is 
consistent with our expectation. There also appears a flat region (minimal-sensitivity 
scale [32|) in the new prediction, which is absent in the former prediction. 

We estimate the error of the new prediction. It is natural to use the scale dependence 
around the minimal-sensitivity scale (~ 160 GeV)Il Following the standard prescription 
we vary the scale by factors 1/2 and 2. When the scale p is varied between 80 and 
320 GeV, the IS energy varies by about 20 MeV below and above the minimal-sensitivity 
scale, respectively. Therefore, the sum of the absolute variations of the IS energy level 
is about 40 MeV|§ The corresponding variation of the top quark MS mass is almost 
one half of it, leading to about 20 MeV, which we take as an error estimate. Another 
error estimate may be obtained from the difference between the NNLO prediction at the 
minimal sensitivity scale (at NNLO) and that at NNNLO, namely the difference between 
the values of Mu{lS) at the local maxima at NNLO and NNNLO in Fig. jS] This gives 
30 MeV as an uncertainty for the top quark mass. For reference, we show the series 
expansion in e at the minimal sensitivity scale at NNNLO: 

Mti{lS) = 2 X (165 + 7.20 + 1.22 + 0.216 + 0.0077) GeV for p = 162 GeV, (7) 

which shows a healthy convergence behavior [mt = 165 GeV and as{Mz) = 0.1185]. 
Thus, we estimate an error in the top quark MS mass determination from Mtt(lS') to be 
20-30 MeV. 

t From the general argument based on the renormalon dominance hypothesis, the minimal-sensitivity 
scale is expected to increase as the perturbation order is raised [25]; see Fig. jS] 

§ This is a factor 2 more conservative estimate than taking the maximal variation of the I/S' energy 
level in this range. 
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We note that the naive error estimate of order [asmY by the uncanceled 

renormalon at u = 3/2 in [30] is order 3-10 MeV, which is still somewhat smaller 
than the current error estimate. This means that the current perturbative order would 
not be high enough to be limited by this renormalon uncertainty. Contribution of the 
renormalon at m = — 1 in the pole mass is estimated to be a few MeV or less (see the 
Appendix), while contribution of the u = +1 renormalon is estimated naively to be order 
Aqcd/’^ ~ 0.5-1.5 MeV (corresponding to Aqcd ~ 300-500 MeV). 

In Ref. [30] the range of the scale variation was taken differently from the above 
range, since no minimal-sensitivity scale for the IS* energy exists for that prediction and 
a different criterion was used. We may check consistency. If we vary the scale in the 
above range for the previous prediction, we obtain the same error estimate for the top 
quark MS mass as in [30] (about 40 MeV). 

Thus, we obtained a better possible accuracy of the top quark mass determination at 
a future linear collider over the previous estimate, which relied only on the renormalon 
dominance hypothesis before the full computations of 03 and ^ 3 . We consider that it 
is not a sheer numerical accident but with a reasoning that we obtain a smaller error 
estimate. Namely, from the general property of QCD a stronger IR cancellation than 
what is predicted by the renormalon dominance hypothesis follows. This interpretation 
is supported by a detailed comparison between the estimates of for rii = 3,4,5 from 
stability of 2 mpoie + Vqcd(?") and the estimates by the renormalon dominance, and also 
by an overall consistent picture drawn in the hrst part of this paper. 

To clarify the current status, we show in Fig. [5] dependences of the IS* energy level 
on the current uncertainty of the exact value of and on the input value of as{Mz) = 
0.1185 ± 0.0006 [33]. The former induces about 10 MeV variation (5 MeV for the top 
quark mass) at the minimal-sensitivity scale, while the latter induces about 90 MeV 
(45 MeV for the top mass) variation. Hence, a precise determination of as(M^), of the 
order of ±0.0001 accuracy, is prerequisite to achieve 20-30 MeV accuracy of the top 
quark mass determination. Prediction of d^ with higher precision is also favorable. 

For comparison, we perform a similar analysis using the potential subtracted (PS) 
mass [3l] as the input parameter. (The dehnition of the NNNLO PS mass is given in [27]. i 
Fig. IHlshows the scale dependence of the toponium IS energy level, where we use the PS 
mass rnps(/i/^ps = 20 GeV) = 173 GeV. To compare with the MS mass, we vary the scale 
from 80 GeV to 320 GeV and hnd the variation of the IS* energy level of about 75 MeV. 
(For 50 GeV < p < 350 GeV, the variation is about 100 MeV, which is consistent with 
[35] .) The uncertainty of as{Mz) causes ± 8 MeV shift of the NNNLO energy level 0 
Thus, use of the PS mass leads to a larger scale variation of the perturbative prediction 
for the IS energy level than the MS mass. We observe qualitatively different scale 
dependences between the two schemes by comparing Figs. [5] and [ 6 ] where this tendency 
is apparent not only at NNNLO but also at lower orders. Furthermore, we conhrm a 
similar tendency in the scale dependences for other n/s, where the values of ds vary 
considerably. We also note that the conversion formula between the PS and MS masses 
induces a scale uncertainty of order 30 MeV for 80 GeV < p < 320 GeV provided that 
mps = 173 GeV is an input value. 

^ The dependence of the PS mass on as{Mz) starts from the order a^, which is the reason for a 
smaller dependence compared to the MS mass. 
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347.5 



Figure 5: Scale dependence of the toponium IS energy level. The input MS mass is taken as 
mt = 165 GeV. Each band for the NNNLO prediction corresponds to variation of inside 

its error (±21.5), where the upper (lower) line in each band corresponds to the upper (lower) 
value of The different bands correspond to different input values of as{Mz)- Predictions 

at lower orders are for as{Mz) = 0.1185. 



H [GeV] 

Figure 6: Scale dependence of the IS energy level in the PS-mass scheme at different orders. 
The input PS mass is taken as mps(/r/,ps = 20GeV) = 173GeV. Two lines for the NNNLO 
result correspond to as{Mz) = 0.1179 and 0.1191. At lower orders as{Mz) = 0.1185 is used. 
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(a) (b) 

Figure 7: Scale dependence of the 15 energy level in the PS-mass scheme at different orders in 
the cases (a) /i/,ps = 50 GeV, mps = 171.2 GeV, and (b) ^/,ps = 80 GeV, mps = 169.5 GeV. 
The values of the PS mass are chosen such that the lines fit in the same range as in Fig. [6l 


Intuitively the difference between using the MS and PS masses may be understood 
as follows. In the MS-mass scheme, the energy of the toponium bound state consists 
of (i) the MS masses of t and i, (ii) contributions to the self-energies of t and t not 
renormalized into the MS mass (typically from gluons whose wavelengths are larger than 
the Compton wavelength of t, Xg ^ and (hi) the potential energy between t and 

t. IR contributions between (ii) and (iii) (typically from Xg larger than the bound-state 
size) get canceled, where the domain of IR cancellation is determined dynamically by 
the wave function of the bound state [9l [36l [26] . 

The composition of the energy of the bound state in the PS-mass scheme is similar, 
except that the renormalized mass (i) is replaced by the PS mass, which renormalizes 
the top quark self-energy from Xg ^ l//i/,ps. In the computation of the self-energy, a 
sharp cut-off is introduced in momentum space at the factorization scale /i/,ps5 which 
is chosen to be of the order of the Bohr scale ~ The cut-off induces a power 

dependence of the PS mass on /i/,ps. Since the l/p/^ps is close to the bound-state size, 
the IR cancellation can become incomplete by artihcial cut-off effects if /i/,ps is too low. 
Such effects tend to be enhanced, due to the increase of the coupling constant at IR and 
the power dependence on /i/,ps- 

We may check consistency of this picture, by conmuting the energy level in the case 
that /U/^ps is taken to be larger than the Bohr scaleHj] In this case, the behavior of the 
predictions in the PS-mass scheme is expected to approach qualitatively that of the 
MS-mass scheme, as only shorter-wavelength contributions are renormalized in the PS 
mass and IR cancellation becomes more complete (artifact of cut-off diminishes). We 
show in Figs. [71(a) (b) the energy level for pyps = 50 and 80 GeV, to be compared with 
Figs. [5l [HI and conhrm this tendency. (We conhrm qualitatively similar behavior for the 
bottomonium energy level as well.) 

Let US discuss other sources of errors. Besides what we have analyzed here, there 

II In principle this is at odds with the standard connting of e in the PS-mass scheme. Furthermore, 
the approximation of subtracting the IR part of the pole mass by an integral of — Fqcd( 9)/2 becomes 
worse as /r/,ps approaches m^. Hence, we take the cut-off in the range asirit < < mt- 
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ni 

0 

1 

2 

3 

4 

5 

6 

I 

II 

31.4 

26.1 

21.3 

17.2 

13.7 

10.6 

8.1 


Table 3: Estimates by the large-/3o approximation for the contribution of the u = —1 UV 
renormalon to ^ 3 . 


are many sources of uncertainties, both of theoretical and experimental origins, in the 
actual top quark mass determination at ILC. Theoretically, these include effects of mixed 
electroweak and QCD corrections (finite width corrections, non-resonant diagrams, non- 
factorizable corrections, etc.), uncertainties in the normalization and shape of the thresh¬ 
old cross section, contributions from higher-spin quarkonium states, method for smooth 
matching to the high-energy cross section, and so forth. In addition effects of the initial- 
state radiation and beam energy spread need to be taken into account in a realistic 
experimental situation for the top quark threshold scan. (See [33 [38] for recent simula¬ 
tion studies for the threshold scan at ILC.) Since feasibility of a high precision top quark 
mass determination can be addressed only by realistic simulation studies incorporating 
all the above effects, the accuracy we present here is what can be achieved in princi¬ 
ple, as a limitation from perturbative QCD. Nevertheless, such a precision is a unique 
possibility achievable only at a future e+e“ collider and worth pursuing. 


Note Added: 

After we completed our work, an analysis was reported on the top quark mass deter¬ 
mination using the NNNLO ti cross section near threshold and using the PS mass [35] - 
Their estimate of about 50 MeV accuracy is larger than the estimate presented in this 
paper (20-30 MeV), which is based only on the uncertainty of the IS* energy level using 
the MS mass. Currently it remains an open question, in the case that the cross section 
is computed thoroughly in terms of the MS mass only, whether the latter estimate is 
increased substantially due to an uncertainty in the shape of the threshold cross section. 
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A UV renormalon 

In this appendix we estimate contributions to the pole-MS mass relation from the UV 
renormalon at u = — 1 using the large-/J q approximation and estimate an uncertainty 
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originating from this renormalon. Using the formula in [18], the contribution to dn 
[defined in eq. ([T])] from the pole at m = — 1 is given by 





(large-/lo approx.), 


( 8 ) 


where Cp = 4/3 is the color factor. In particular the contributions to are evaluated 
explicitly for various rii in Tab. [3l Comparing these values with the corresponding errors 
of in Tab. O we hnd that they are smaller than the errors by factors 4-6 for 

0 < rii < 4 and by factors 10-20 for ni = 5,6. This is consistent, since each error of 
dfm is dominated by the contribution from the analytic part at m = 1/2 and 
belongs to the analytic part. It suggests that the u = —1 UV renormalon is not a 
dominant component of the contribution from the analytic part (for n = 3). 

In the rest of this appendix we estimate the contribution of the u = —1 UV renor¬ 
malon to the quarkonium IS* energy level, taking the toponium case (n/ = 5) as an 
example. (The case for the bottomonium is qualitatively similar.) 

In the IS* energy level (at the leading-logarithms) only the pole mass contains the 
u = —1 UV renormalon. In general a UV renormalon induces a factorial growth of 
perturbative series, as shown in eq. dHj) (similarly to an IR renormalon), which breaks 
convergence of the perturbative series. Nevertheless, since the corresponding singularity 
in the Borel plane (w-plane) lies along the negative real axis, a dehnite value can be 
assigned to the contributions of a UV renormalon by Borel summation. The perturbative 
series corresponding to a UV renormalon converges up to a certain order (n < n*) and 
diverges beyond that order {n > n*), which is a typical feature of an asymptotic series. 
In the case of the u = —1 UV renormalon (the UV renormalon nearest to the origin in 
the Borel plane), the critical order n* is given by 


An 

^oas{mt) 


^ 15. 


(9) 


Therefore, the perturbative series is still converging in our NNNLO calculation. The 
hrst several terms of the u = —1 contribution (in the large-/S q approximation) read 


= 2mt 


1 


u=-l] 


n=0 


a 


_ \ \ n+l 




( 10 ) 


= 2 X [165 - 1.44 + 0.096 - 0.013 + 0.0025 - 0.00068 + •••], (11) 


for jj, = rut = 165 GeV and as{mt) = 0.109. According to a standard estimate with 
an asymptotic series, the error of the prediction is of the order of the last known term. 
Hence, at NNNLO, we can estimate the error due to the u = —1 renormalon to be of 
order 2.5 MeV for the top quark mass determination. 

Alternatively we can estimate the error using the difference between the Borel summed 
value and the perturbative contribution up to NNNLO: 


6mt 


ACpfut r, 

/ du 

esfjf) Jo 


1 

1 + u 


{1 — u + — u^) 


exp 


Anu 

(Joasimt) 


^ —0.51 MeV. 


( 12 ) 
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It is somewhat smaller than the above estimate. [An error estimate by the N^LO term 
of eq. (ITT]) gives a better estimate.] 

From the above examinations, one expects that the contribntion of the u = —1 
renormalon is fairly modest and minor in the error estimate in the determination of the 
top qnark mass, which is performed in the main body of this paper. As long as the 
perturbative series is converging, the error due to the u = —1 renormalon decreases. 
This is in contrast to the n = 3/2 renormalon, which induces a limitation in achievable 
accuracy of order ■ A crude estimate based on the large-/9o approximation 

indicates that at NNNLO the error due to the u = —1 renormalon is smaller than the 
error due to the u = 3/2 renormalon. 
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